The dimer-hole-RVB state of the 2-leg t-J ladder: A Recurrent 

Variational Ansatz 

German Sierra* 

Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106 

Miguel Angel Martm-Delgado^ 
Departamento de Fisica Teorica I, Universidad Complutense, 28040 Madrid, Spain 

Jorge Dukelsky-t 

Instituto de Estructura de la Materia, C.S.I.C., 28006 Madrid, Spain 

Steven R. White^ 

Department of Physics and Astronomy, University of California, Irvine, CA 92697 

D.J. Scalapino** 

Department of Physics, University of California, Santa Barbara, CA 93106 

(February 7, 2008) 

Abstract 

We present a variational treatment of the ground state of the 2-leg t-J ladder, 

*0n leave from Instituto de Matematicas y Fi'sica Fundamental, C.S.I.C., 28006 Madrid, Spain, 
em: sierra@sisifo.imaff.csic.es 
^em: mardel@eucmax.sim.ucm.es 
•fem: emduke@iem.csic.es 
§em: srwhite@uci.edu 
**em: djs@spock.physics.ucsb.edu 

1 



which combines the dimer and the hard-core boson models into one effective 
model. This model allows us to study the local structure of the hole pairs as a 
function of doping. A second order recursion relation is used to generate the 
variational wave function, which substantially simplifies the computations. 
We obtain good agreement with numerical density matrix renormalization 
group results for the ground state energy in the strong coupling regime. We 
find that the local structure of the pairs depends upon whether the ladder is 
slightly or strongly dopped. 
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Typeset using REVT^ 



INTRODUCTION 



The 2-leg, t — J ladder represents one of the simplest systems which exhibits some of 
the phenomena associated with high Tc cuprate superconductivity The ground state 

of the undoped system, a 2-leg Heisenberg ladder, is a spin liquid with a finite spin gap and 
exponentially decaying antiferromagnetic spin-spin correlations. Upon doping, the spin gap 
remains and there appear power law CDW and singlet superconducting pairing correlations. 
In addition, the pairing correlations have an internal d^^i-y^-Wke symmetry with a relative 
sign difference between the leg and rung singlets which make up a pair. Despite all of the 
numerical and analytical work which has been done on this system, we still lack a picture of 
the ground state which accommodates all of these physical properties. There are, however, 
many hints of what that picture may look like. It is the purpose of this paper to take one 
step further in that direction. 

Short-range resonating valence bonds (RVB) provide a useful basis for representing the 
ground state of spin liquids [^If. For the t — J ladder, a O^'^-order picture has been provided 
by the study of the strong coupling limit where the exchange coupling constant along the 
rungs, J', is much larger than any other scale in the problem. The other coupling constants 
of the model are, J: the exchange coupling constant along the legs, and t and t': the hopping 
parameters along the legs and the rungs respectively. In the limit J' » J,t,t', the ground 
state of the undoped ladder is simply given by the coherent superposition of singlets across 
the rungs. Addition of one hole requires the breaking of one of these singlets, in which case 
the hole gets effectively bound to the unpaired spin, becoming a quasiparticle with spin 1/2 
and charge |e|. Addition of another hole leads to the binding of two holes in the same rung 
in order to minimize the cost in energy. In this picture there is no spin-charge separation, a 
fact that remains valid down to intermediate and weak couplings, as confirmed by various 
numerical and analytical studies. Based on this picture it is possible to construct an effective 
theory describing the motion and interactions of the hole pairs . It is given by a hard-core 
boson model (HCB) characterized by an effective hopping parameter t* and interaction V* 
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of the hole pairs. The HCB model describes the doped ladder as a Luther-Emery liquid, 
with gapped spin excitations and gapless charge collective modes, which are responsible for 
the CDW and SC power law correlations. We summarize the O^'^-order picture in Figure 1, 
which shows a typical state of HCB's, as well as the two building blocks that are used its 
construction. 

In order to go beyond this picture, we need to consider the fluctuations of the states 
of the HCB model. To lowest order in perturbation theory they are shown in Fig. 2. The 
admixture of the state shown in Fig. 2(a) is of order J / J' and represents a resonance of 
two nearest neighbor rung singlets. According to the standard RVB scenario, this resonance 
effect leads to a substantial lowering of the ground state energy. The state in Fig. 2(b) 
is of order t/J', and it can be though of as a bound state of two quasiparticles, whose 
characteristic feature is the diagonal frustrating bond across the holes. From the RVB point 
of view. Fig. 2(b) is a resonance of a singlet and a hole pair. The importance of this state, 
even for intermediate couplings such as J = J' = 0.5t, was emphasized in the DMRG study 
of reference , where it was shown to be the most probable configuration of two dynamical 
holes in a 2- leg ladder. In the HCB model of ||^, the states of the form of Fig. 2(b) are 
taken into account as intermediate or virtual states, which lead to the effective hoping, t* 
and interaction, V* between the hole pairs. It is clear however that "integrating out" the 
diagonal states through perturbation theory, erases the internal structure of the hole pairs. 
Here we want to extend the HCB description to include the internal structure of the hole 
pairs. 

In order to define an effective model which would retain the degrees of freedom associated 
with the internal structure of the hole pairs, we need to consider the states that appear in 
second order in the strong coupling expansion. They are given in Fig. 3. Let us comment 
on them. The state of Fig. 3(a) is of order (J/ J'Y and it is a higher order RVB state, whose 



contribution to the ground state of the undoped ladder was studied in . In this reference 
it was shown that its inclusion in a variational ansatz improves the numerical results, but 
does not change the qualitative picture obtained using the dimer ansatz BJTl||. The state of 



Fig. 3(b), which is in fact first order in t', can be seen as a bound state of two quasiparticles, 
while 3(c) and 3(d) are higher order corrections to the diagonal state shown in Fig. 2(b). For 
these reasons it seems consistent to keep the state 3(b) on an equal footing with the states 
2(a) and 2(b). To give further support to this choice, we notice that the exact solution for 
two holes on the 2x2 cluster requires a superposition of the states shown in 2(b) and 3(b) 
along with 1(a) and 1(b) (see Fig. 4) 

In summary, we conjecture that in order to discuss the nature of the superconducting 
order parameter of the doped 2-leg, t — J ladder, in the strong coupling regime, it is sufficient 
to consider states built up from 5 possible local configurations, given by rung-singlet-bonds 
(Fig. 1(a)), rung-hole-pairs (Fig. 1(b)), two-leg-bonds (Fig. 2(a)), hole-pairs with a singlet 
diagonal bond (Fig. 2(b)) and hole-pairs with a singlet leg bond (Fig. 3(b)). A typical state 
constructed using these building blocks is shown in Fig. 5. We shall call these types of states 
dimer-hole-RVB states. The effective model that governs their dynamics will be called the 
dimer hard-core boson model (DHCB) and its Hamiltonian can be determined by considering 
the fiuctuations of the dimer-hole states, in a manner similar to the one considered above for 
the HCB states. The DHCB model contains spin and charge degrees of freedom, together 
with their couplings, and in that sense is an interesting model to study the interplay between 
the two types of degrees of freedom, although here we will focus on the variational ground 
state of the model. 

The mathematical formulation of the DHCB model involves an interesting but com- 
plicated combination of vertex and Interaction Round a Face (IRF) models. The latter 



terminology is borrowed from Statistical Mechanics [|T^]. The vertex variables describe the 
number of electrons per rung, i.e. = 0, 1, 2, while the IRF variables describe the number 
and type of bonds connecting two rungs, i.e. ii^i+i = 0, 1^, Ih, 2, where the subindices d, h 
indicate the diagonal or horizontal nature of the bond. The only allowed configurations for 
two consecutive IRF variables (£i,i+i, £1+1,4+2) are: (0, 0), (1^, 0), (1/j, 0), (2, 0) together with 
their permutations. Moreover the vertex variables are subject to certain constraints imposed 
by the IRF ones. Namely, A) if £i,i+i = Id or 1^ then Ui = nj+i = 1, and B) if £j,j+i = 2 



then Hi = rii+i = 2. Only if = can rii and n^+i take any value, i.e. 0, 1 or 2. 

It is beyond the scope of this work to present a full account of the DHCB model. Instead, 
we shall try to uncover some of its physics, by means of a combination of two approaches, 
namely the Density Matrix Renormalization Group |13[ and the Recurrence Relation Method 
(RRM) |T^. While the DMRG is a powerful numerical technique, which in many cases 
yields the exact answer, the RRM is essentially analytic, lacking the numerical precision 
of the DMRG, but sharing with it some features, as for example the Wilsonian way of 
growing the system by the addition of sites at the boundary. In the RRM one begins with 
an assumption about the local configurations through which the system grows. Then one 
may test whether the state that is generated gives results in agreement with the essentially 
exact DMRG results. 



THE VARIATIONAL WAVE FUNCTION 

The Hamiltonian of the 2-leg, t — J ladder is given by, 

H = Hs + Hk = Jij (Sj ■ Sj — \ninj) 

- Y.{i,j),sUj Pg {c\,sCj,s + c]^gCi^s) Pg (1) 

where J^j, tij = J,t or J', t', depending on whether the hnk (ij) is along the legs or the rungs 
respectively. Pq is the Gutzwiller projection operator which forbids double occupancy. The 
rest of the operators appearing in (P are standard (we use the conventions of reference 0). 
Each site i is labelled by the coordinates (x, y) with x = 1, . . . , N and y = 1,2. We choose 
open boundary conditions along the legs of the ladder. 

The pair field operator which creates a pair of electrons, at the sites i and j, out of the 
vacuum is given by, 

^. = ^i<^,Ai + 4Ai) (2) 

As explained in the introduction, we want to built up an ansatz for the ground state 
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based on the 5 local configurations of the DHCB model. The explicit realization of these 
configurations in terms of pair field operators are given by ( see Fig. 6), 

\Mx= |o). 

|01,o)x = ^(x,l){:.,2) |0)^ 
\<p2fi)x,x+l = -U ^\x,2)(x+l,2) 10)^.2^+1 

\<p2,l)x,x+l = [b i^lx,l)(x+l,2) + ^(x,2)(z+l,l)) 

+ ^(a;,2)(x+l,2))] \^)x,x+l 

where \0)x is the Fock vacuum associated with the rung labelled by the coordinate x ( 
\0)x,x+i = |0)x ® |0)x+i)- The states |0„,p), involve n = 1,2 rungs and p = 0,1 pairs of 
holes. The variational parameter u gives the amplitude of the resonance of a pair of bonds 
between vertical and horizontal positions while b and c are the variational parameters 
associated with the diagonal and horizontal configurations of two holes respectively. In the 
strong coupling limit, J' » J,t,t' , we expect to find u J / J' , 6 ~ tj J' and c ~ tt' / J'"^. 

Let us call \N,P) the ground state of a ladder with rungs and P pairs of holes. Of 
course we should be in a regime of the coupling constants where there is binding of two 
holes. The state \N,P) will be in general a linear superposition of the dimer-hole states 
of Fig. 5, which suggests that working with this sort of states could be a formidable task. 



Fortunately, we can apply the method developed in [|T0[ to generate \N,P) in a recursive 
manner, in terms of the states of the ladders with — 1 and N — 2 rungs, and P and P — 1 
pairs of holes. In it was shown that |A^, P = 0), which is in fact a dimer-RVB state 0,|11 



can be generated by a second order recursion relation. Then by a simple procedure one can 
compute overlaps and expectation values of different operators using recursion formulas, 
whose thermodynamic limit can be studied analytically. 

Following the strategy of considering first the HCB states and then the DHCB ones, we 
shall give the rule that generates the former type of states. It is given by the first order 
recursion relation, 

\N + 1,P + 1) = \N,P + 1) |0i,o)iv+i + \N, P) \<Pi,i)n+i (4) 



supplemented with the initial conditions, 

|l,0) = |</)i,o) 

|l,l) = l0i,i) (5) 
\N,P) = 0, hi N < P 

Calling F^p^ the number of linearly independent states contained in \N,P), we deduce 
from Eq.(|^) the recursion relation, 

T^HCB _ T^HCB I r^HCB ff^\ 
^N+1,P+1 — ^N,P+1 "T ^N,P \^) 



whose solution is given by the combinatorial number, 

(7) 



HCB 



TP 



Eq. is the dimension of the Hilbert space of the HCB model with sites and P pair 
of holes. We have not introduced variational parameters in Eqs. (|^), but if we did, then all 
states of the Hilbert space of the HCB model would be generated by the first order recursion 
relation. It may be worthwhile to recall that the HCB model is essentially equivalent to the 
spinless fermion model or the XXZ model 0. 

Turning now to the DHCB model, the key point is to realize that the dimer-hole states 
can be generated by the following second order recursion relation, involving the local con- 
figurations given by eq.(||), 

|iV + 2, P + 1) = |iV + 1, P + 1) |0i,o)^+2 

+ \N + 1,P) 

N+2 + \N,P+1) 102 ,q) N+l,N+2 + \N, P) \4'2,i)n+1,N+2 (8) 

with the initial conditions (^). See Fig. 7 for a graphical representation of (P). 
Counting dimer-hole states 

Let P/v,p denote the number of dimer-hole states of a 2-leg ladder with N rungs containing 
P pairs of holes. According to (|^) they satisfy the recursion relation 
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Fn+2,P+1 

with the initial conditions 



— -^Af+i,p+i + Fn,p+i + Fjsf^i^p + AFn^p 



(9) 



Fn,n = 1, Fn,p = for N <P (10) 
From (PI) and (pOl) we deduce that F/vo satisfies the well known Fibonacci recursion 



formula |10|, and that in the limit of very large N it grows exponentially, 



'N,0 



(iV»l) (11) 



where $o = |(1 + V^) is the golden ratio. Using generating function methods |jTO[ one can 



easily solved the recursion relation (|^), together with the initial condition (p!oD. The result 
is given by the contour integral, 

dz z^+i(2 + 4)^ 



(12) 



2m {z^-z- 1)^+1 

where the contour encircles the singularities of the integrand. For P = the integrand has 
two simple poles at the zeros of the polynomial — z — 1, the largest of which is precisely 
the golden ratio $o- In this way one gets Eq.(0). For a finite number of holes the residue 
formula applied to (|l^) yields, to leading order in N 

Fn,p ~ A^-^ N»l, P : finite (13) 

where the proportionality constant depends only on P. Let us finally consider the limit 
where both and P go to infinity, while keeping their ratio fixed. 

Number of holes -P „ . 

^ = T-T — r-- — = T7 > < X < 1 (14) 

Number of sites iV 

Here x can be identified with the hole doping factor of the state \N,P). The saddle point 
method applied to ([T^) gives the asymptotic behaviour of the number of dimer-hole states 
for a finite density of holes, 

F.,.~f(.r, /w = ^^ii±iii, (15) 
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where $ = is the highest root of the following equation 

^ -$-!)($ + 4) 

$($2 + 8$ -3) ^ ^ 

The function /(x) is depicted in Fig. 8. Observe that $(0) = $o- The effect of a 
finite density of holes is that of moving a singularity. This phenomena also occurs in the 
computation of the energy, and other observables. 

Ground State Energy 

The parameters u, b, c are found by the standard minimization of the mean value of the 
energy (A^, P\Hn\N, P) / {N, P\N,P), where if at denotes the Hamiltonian of the ladder with 

rungs. The usefulness of Eq.@ is that it implies that the wave function and energy 
overlaps also satisfy recursion relations. Let us define the following quantities, 

Z^,P= {N,P\N,P) 

Yn,p= n{<Pi,o\{N-1,P\N,P) 

Ej,^P= {N,P\H^\N,P) (17) 
Dn,p= N{h,o\{N-l,P\H^\N,P) 
Wm,p= {N,P\nM\N,P) 

where is the number operator acting on the rung A^. The off-diagonal overlaps arise from 
the cross terms when applying (§) to the ket and the bras in (A^ + 2, P + 1|A^ + 2, P + 1) 
and (A^ + 2, P + l\Hi\f^2\N + 2, P + 1). The recursion relations satisfied by (|T7D are given 

by, 
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Zn+2,P+1 = -^AT+l.P+l + Zn,P+1 + U Fat+i^p+i + Z^^i^p + 2(6^ + C^) Zn^p 

-E'Ar+2,P+l = -E'a^+i.p+i — J' Zn+1,P+1 + Em,p+i — (2 J + J' /2)u'^ Zm,P+1 + En+1,P 
+2(62 + c2)E7v,p - (2 Jc^ + 46t + 8bct') Zn,p + uDn+i,p+i 

-2u{J + JV2)FAr+i,p+i - 4t6FAr+i,p (18) 
W^7V+1,P+1 - ijMW;v,P+l - + c'')Wm,P 

Dn+2,P+1 = -E'a'+I.P+I — J' Zn+1,P+1 + w/2DAr+i_p+i — u{J + J'/2) Y/v+l,P+l 

—2tbZN^P — ^JWn+i^p^i 
^Af+2,p+i = 2Z7v+i,p+i + 2^2 Zjv^p+1 + 2(6^ + c^)Z]\f p + 2My)v+i^p+i 

The initial conditions read, 

Zo,o = 1, ^0,0 = -£'0,0 = Dqo = Wo ft = 

(19) 

Xjv,p = 0, for N <P and X = D, W 

For finite values of N and P, and given choices of u, b, c, one can iterate numerically 
the recursion relation using the initial conditions ( ]TP| ) and look for the minimum of the 
ground state energy Ejq^p / Z^^p. We give below the results obtained using this variational 
method for a 2 x 32 ladder and compare them with the corresponding results obtained with 
the DMRG. 

THE RRM WAVE FUNCTION VERSUS THE DMRG: NUMERICAL RESULTS 

As explained in the introduction the DHCB model is the appropiate framework to study 
the strong coupling limit of the 2-leg ladder, if one wishes to take into account the local 
structure of the hole pairs. To check the validity of this asumption we have studied the 
cases where the coupling constants takes the following values, t = t' = 1, J = 0.5 and 
J' = 0.5, 1,2,3,4 and 5. In this manner we go from the intermediate coupling regime, i.e. 
J' ~ 1 to the strong coupling regime J' ^ 3.5. We are always working in a non-phase- 
separated region. 
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In Figure 9 we show the ground state energy of the 2 x 32 ladder, for the previous choices 
of parameters, computed with the RRM for all dopings and the DMRG for x = 1/8, 1/2 and 
7/8. One sees that the results obtained with the RRM wave function agree reasonably well 
with those of the DMRG and their accuracy improves as J' increases. 

The kinetic energy of the ladder is shown in Fig. 10. It has the pattern expected for a 
collective charge mode, as described by the HCB and the DHCB models. The similarity 
between this figure and Fig. 8 have a common origin. They both correspond to holes moving 
collectively through the spins in a complicated many body state. 

Fig. 10 shows the existence of an optimal doping for which the kinetic energy is a min- 
imum. The existence and position of this optimal doping depends on the values of the 
coupling constants. 

The nature of this many body state is clarified by figures 11, 12 and 13 where we show 
the values of the variational parameters u, b and c as functions of the doping x for different 
coupling constants. The parameter u starts from a positive value corresponding to the 
undoped ladder and it decreases upon doping until a critical value Xc{J/J'), where it 
vanishes. For higher dopings u becomes negative. For the undoped ladder the parameter u 
can be interpreted as the square of the RVB amplitude /iRVB for having a bond along the 
legs [l^. The analogue amplitude for a bond along the rungs has been implicitly normalized 



to 1. For low doping, i.e. x < Xc, since u{x) > 0, we can similarly define a doping dependent 
amplitude for a leg-bond as 

u{x) = /irvb(^) > 0> < ^c) (20) 

In order to fulfill the Marshall theorem for the undoped ladder one requires the RVB am- 
plitude /irvb(O) to be positive |Q, which explains why u(0) is also positive. Actually for the 
positivity of u(0) one just need /irvb(O) to be a real number. At x = /iRVB(O) increases with 
J/ J' due to the resonance between rung and leg singlets, according to the RVB scenario. 
Upon doping, however, the holes give rise to destructive interference which degrades pro- 
gressively the aforementioned resonance mechanism. This explains why u{x) and /iRVB(a^) 
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decrease with x. For x < Xc the ground state is dominated by the resonating valence bonds 
and the RVB picture remains quahtatively correct. 

For X > Xc the interference due to the holes has driven u negative and it is no longer 
appropiate to interpret u{x) as the square of /irvb- Rather, the physical interpretation of 
the overdoped region comes from the solution of the Cooper problem in the t — J, 2-leg 
ladder, and its BCS extension. It can be shown analytically that two electrons in the latter 
system form a bound state only under certain conditions (details will be given elsewhere). 
For J = 0.5, t = t' = 1 one must have J' > 3.3048, ( note that the binding of two electrons in 
the t — J chain requires J/2t > 1 [Q)- The exact solution for 4 or more electrons is difficult 
to construct, but we expect it to be given essentially by a Gutzwiller projected BCS hke 
wave function. A short range version of the latter type of wave function can be generated 
from the recursion relation (H), with u a negative parameter, which can be written as 

u{x) = -hlcsix) < 0, (x > Xc) (21) 

where h^cs is the BCS amplitude for finding two electrons at distance 1 along the legs. 
Of course this interpretation of u as minus the square of a BCS amplitude requires it to 
be negative. As we put more electrons into the ladder the value of /iscs decreases and for 
electron densities larger than 1 — Xc, we switch into the RVB regime. 

The difference between the underdoped and overdoped regimens can be attributed to 
two different internal structures of the pairs. In the low doping regime x < Xc, holes doped 
into the spin-liquid RVB state form pairs with an internal (i^2_y2-like structure relative to 
the undoped system. However for x > Xc one moves into the low density limit characterized 
by electrons doped into an internal s-wave like symmetry. This issue will be discussed in 
detail in a separate publication. 

Let us now comment on Figs. 12 and 13. Both are very similar and show that for x ~ 1/2, 
b and c reach their maximum. At x = 1/2 there are as many electrons as holes, and in a 
certain sense the ground state of the ladder is a large scale reproduction of the microscopic 
ground state of the 2x2 cluster given in Fig. 4. Indeed for J = J' = 0.5, t = t' = 1 the ratio 
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h/a oi the parameters appearing in Fig. 4 is given by 1.30, which is very close to the value of 
h at its maximum. For x < 0.7 and J' = 0.5 the parameter b is larger than 1 and it is always 
larger than c for all dopings and couplings. This is in agreement with the DMRG results 
of 1^, which show the importance of the diagonal frustrating bonds above the horizontal or 
vertical ones for J/t = J' /t = 0.5. 

Finally Fig. 14 is a J/t — n diagram which shows the boundary of phase separation 
obtained by means of the DMRG and the RRM in the case where J = J' ,t = t' = 1. 
Observe that this is not the strong coupling case we have been discussing so far, and hence 
the validity of the RRM is more questionable. In any case, we see an overall agreement 
between both results (see references p,|15|,|16[ for comparisons with other numerical results). 
In the two-leg t-J model, phase separation is controlled by J, rather than J', so the strongest 
coupling we have considered above, J'/t = 5, J/t = 0.5, t'/t = 1, does not phase separate. 



CONCLUSIONS 

In this paper we have proposed an extension of the effective hard-core boson model 
(HCB) of the 2-leg ladder of reference P], in order to include the local structure of the 
hole pairs. The extended effective model, called the DHCB model, contains both dimer 
bonds, hard core bosons and various combinations between bonds and holes, whose relevance 



have been studied previously with DMRG p. Generalizing the methods of reference IjlO 
to the case with holes, we study a variational ansatz for the ground state of the DHCB 
model, which depends only on three variational parameters. The resulting dimer-hole state 
is generated by a second order recursion formula, which also leads to recursion formulas 
for the overlaps necessary to compute the energy of the ansatz. We give the results of the 
energy minimization for the 2 x 32 ladder and compare them with those obtained with 
the DMRG method in the strong coupling region. The recursion relations we have derived 
for the ground state energy can be solved analytically in the thermodynamic limit and the 
minimization can be then done numerically. Finally we give a physical interpretation of the 
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behaviour of the variational parameters with doping. 



ACKNOWLEDGMENTS 

GS would like to thanks the organizers of the ITP program "Quantum Field Theory in 
Low Dimensions: Prom Condensed Matter to Particle Physics" for the warm hospitality. 
MAMD thanks the organizers of the Benasque Center of Physics 1997 for their support 
and hospitality. GS acknowledges support from the NSF under Grant No. PHY94-07194 
and the Direccion General de Ensenanza Superior, MAMD acknowledges support from the 
CICYT under contract AEN93-0776, JD acknowledges support from the DIGICYT under 
contract No. PB95/0123, SRW acknowledges support from the NSF under Grant No. DMR- 
9509945, and DJS acknowledges support from the NSF under Grant numbers PHY-9407194 
and DMR-9527304. 



15 



REFERENCES 

[1] E. Dagotto, J. Riera and D. J. Scalapino, Phys. Rev. B 45, 5744 (1992). 

[2] T. Barnes, E. Dagotto, J. Riera and E. S. Swanson, Phys. Rev. B 47, 3196 (1993). 

[3] S. Gopalan, T.M. Rice and M. Sigrist, Phys. Rev. B 49, 8901 (1994). 

[4] S.R. White, R.M. Noack and D.J. Scalapino, Phys. Rev. Lett. 73, 886 (1994). 

[5] H. Tsunetsugu, M. Troyer, T. M. Rice, Phys. Rev. B 49, 16078 (1994). 

[6] M. Troyer, H. Tsunetsugu and T. M. Rice, Phys. Rev. B 53, 251 (1996). 

[7] S. A. Kivelson, D.S. Rokhsar and J.P. Sethna, Phys. Rev. B 35, 8865 (1987). 

[8] S. Liang, B. Doucot and P.W. Anderson, Phys. Rev. Lett. 61, 365 ( 1988). 

[9] S.R. White and D.J. Scalapino, Phys. Rev. B 55, 6504 (1997). 

[10] G. Sierra and M.A. Martfn-Delgado, "The Short Range RVB State of the Even Spin 
Ladders: A Recurrent Variational Approach", |cond-mat / 9704212| , Phys. Rev. B , to 
appear. 

[11] Y. Fan and M. Ma, Phys. Rev. B 37, 1820 (1988). 

[12] R. Baxter, "Exactly Solved Models in Statistical Mechanics" , Academic Press, London, 
(1982). 

[13] S.R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. B 48, 10345 ( 1993). 

[14] M. Ogata, M.U. Luchini, S. Sorella and P.P. Assaad, Phys. Rev. Lett. 66, 2388 (1991). 

[15] C.A. Hayward and D. Poilblanc, Phys. Rev. B 53, 11721 (1996). 

[16] K. Sano, Jour. Phys. Soc. Jap. 65, 1146 (1996). 



16 



FIGURE CAPTIONS 



Figure l:The 0**^ order picture of tlie Hard Core Boson model: a) The vertical bond, b) 
the vertical hole-pair singlet, c) a typical state of the HCB model. 

Figure 2 : The two lowest order states in the strong coupling limit J' ^ J, t, t' of 
the HCB model. They represent the first order contribution to the DHCB model, a) the 
resonance of two vertical bonds, b) bound state of two quasiparticles. 

Figure 3 : Higher order strong coupling states contributing to the DHCB model, a) 
a higher order RVB state, b) a bound state of two quasiparticles, c) and d) higher order 
corrections to the diagonal state 2(b)). 

Figures 4 : The exact ground state for a single plaquette with two holes |^ (case N = 2 
and P = 1). 

Figure 5 : A typical dimer- hole- RVB state. 

Figure 6 :Elementary building block states of the RRM used in the construction of the 
dimer-hole states. 

Figure 7 :A pictorical representation of Eq. (^. 

Figure 8 :The function f{x) appearing in (plSj). The maximum appears at x = 0.44. 

Figure 9 : Ground state energy per site of the 2 x 32 ladder with J = 0.5, t = t' = 1 
and J' = 0.5, 1, 2, 3, 4, 5. The remaining data given below in figures 10-13 also corresponds 
to these choices of couplings. The continuum curves are obtained with the RRM, while the 
special symbols are the DMRG data corresponding to x = 1/8, 1/2 and 7/8 respectively. 

Figure 10 : Kinetic energy per site. 

Figure 11 : The variational parameter u as function of the doping. 
Figure 12 : The variational parameter b as function of the doping. 
Figure 13 : The variational parameter c as function of the doping. 
Figure 14 : Boundary of the phase separation region in the case where J = J' ,t = t', 
computed with DMRG and the RRM. 
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